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ABSTRACT 

Surface impedance boundary conditions are employed to reduce 
the solution volume during the analysis of scattering from lossy 
dielectric objects. In a finite difference solution, they also can 
be utilized to avoid using small cells, made necessary by shorter 
wavelengths in conducting media throughout the solution volume. 
The standard approach is to approximate the surface impedance over 
a very small bandwidth by its value at the center frequency, and 
then use that result in the boundary condition. In this paper, two 
implementations of the surface impedance boundary condition are 
presented. One implementation is a constant surface impedance 
boundary condition and the other is a dispersive surface impedance 
boundary condition that is applicable over a very large frequency 
bandwidth and over a large range of conductivities. Frequency 
domain results are presented in one dimension for two conductivity 
values and are compared with exact results. Scattering width 
results from an infinite square cylinder are presented as a two 
dimensional demonstration. Extensions to three dimensions should 
be straightforward. 

I. Introduction 

The Finite Difference Time Domain ( FDTD ) technique permits the 
analysis of interactions of electromagnetic waves with objects of 
arbitrary shape and material composition. This method was first 
proposed by Yee [1] for isotropic, non— dispersive materials in 
1966; and through various modifications during the past twenty 
years, it has evolved into a mature computational technique. 
Reference [2] and the references contained therein provide an 
account of various extensions and modifications of the original 
FDTD algorithm. The present FDTD technique is capable of 
electromagnetic scattering analysis from objects of arbitrary and 
complicated geometrical shape and material composition over a large 
band of frequencies. This technique has recently been extended to 
include dispersive dielectric materials [3], chiral materials [4] 
and plasmas [5]. Due to these numerous capabilities, the FDTD 
method has begun to gain widespread acceptance as a viable 
computational alternative to the classical method of moments (MM) 
technique for many problems. 

To analyze electromagnetic field interaction with lossy 
dielectric objects, the FDTD method requires that the interior of 
the object be modeled in order for fields to penetrate the body. 
Accurate modeling often requires a very fine spatial grid resulting 
in a relatively large number of cells for moderately sized objects. 
A highly conducting dielectric object can be replaced by a surface 
impedance boundary condition (SIBC) that is a function of the 
material parameters. Thus, this boundary condition eliminates the 
spatial quantization of the object and reduces the overall size of 
the solution space not only by eliminating cells within the lossy 
dielectric, but also by allowing larger cells to be used in the 
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exterior region. As with any computational electromagnetic tool, a 
technique that reduces the solution space or number of unknowns is 
quite welcome. 

Of historical interest, surface impedance boundary conditions 
were first proposed by Leontovich in the 1940's [6] and were 
rigorously developed by Senior in 1960 [7]. During the past thirty 
years, researchers have applied surface impedance concepts in the 
frequency domain to numerous electromagnetic scattering problems. 
Time domain surface impedance concepts received little attention 
until recently. Through some impressive work, Maloney and Smith 
[8] have previously implemented a surface impedance boundary 
condition in the FDTD method. However, their implementation has a 
minor disadvantage because the exponential rates and coefficients 
for recursive updating have to be reevaluated each time the 
conductivity or loss tangent is changed. With our proposed method, 
the exponential rates and coefficients only have to be evaluated 
once . Tesche [9] has also investigated surface impedance concepts 
in an integral equation time domain solution, but presented limited 
computational time domain results. 

It is the purpose of this paper to introduce a constant 
surface impedance boundary condition that is applicable for a 
single frequency and a dispersive surface impedance boundary 
condition that is applicable over a large frequency bandwidth and 
range of conductivities. The dispersive surface impedance includes 
frequency variations which results in a time domain boundary 
condition involving a convolution. We will then show how to 
efficiently evaluate this convolved surface impedance using 
recursion. 

II. Motivation 

The motivation for implementing a SIBC in the FDTD method is 
to reduce the computational resource requirements for modeling 
highly conducting lossy dielectric objects. In the standard FDTD 
method, modeling highly conducting lossy dielectric objects 
requires that the cell size be chosen small enough to resolve the 
field inside the object at the maximum frequency of interest. For 
example, suppose scattering from a lossy dielectric object with 
permeability /x, permittivity e and conductivity cr=2 . 0 S/m is to be 
studied over the frequency band 0—10 GHz. The cell size must be 
chosen as some fraction of the wavelength inside the conducting 
material at the maximum frequency of interest. Thus the cell size 
is chosen (typically) as 



Sx = Sy = Sz = — 
10 


( 1 ) 
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where £ r is the complex relative permittivity constant of the 

material and A. and A 0 are the wavelengths inside the material and 
in free space at 10 GHz, respectively. The complex permittivity 
for lossy dielectrics in the frequency domain is 

£ = e + — (2) 

jo 


where o> is the radian frequency. The complex relative permittivity 
is determined using (2) as 



€ r + 


a 

joe 0 


(3) . 


If the material is a good conductor over all frequencies of 
interest, then the constitutive parameters satisfy the condition 

— » 1 (4) . 

(06 


Therefore, £ r can be approximated as 


a * _JL_ 


(5) . 


Assuming parameters m=Mo and e=e 0 and using the values of £ r and A. 0 

at 10 GHz, the cell size is Sx = Sy = Sz = 1.582 mm. If a SIBC is 
used, then the cell size need only be chosen to resolve the field 
in free space and (1) is modified to 

^ x sibc = ^Ysibc = ^ z sibc " (6) 


Again, using the value for A 0 at 10 GHz, the cell size is £x S]BC - 
^Vsibc = 5z sibc = 3 -° nun * Thus the cel1 size has been increased b y 
the factor JYK~\ =1.90, and the computational storage requirements 

are reduced by the same factor. Therefore, the computational 
savings, denoted by S, is 


S 



(V) 


where e r is given by (5), (4) is satisfied for all frequencies of 

interest and d is the number of dimensions. 
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III. FDTD Constant Surface Impedance Implementation 

To implement the constant SIBC in the FDTD method we consider 
the planar air-lossy dielectric interface as shown in Figure 1. 
The conducting material has permittivity e, permeability ji and 
conductivity a. We assume that the thickness of the material is 
large compared to the skin depth. We will also assume that the 
material is linear and isotropic and a basic familiarity with the 
Yee algorithm [1]. Figure 1 also shows the one-dimensional FDTD 
grid. 

The first order (or Leontovich) impedance boundary condition 
relates tangential total field components and is given in the 
frequency domain as [6] 

E x (u) = Z s (w)H y (co) (8) 


where z s (w) is the surface impedance of the conductor. The 
frequency domain surface impedance for good conductors is 


Z s (6>) 


d + j) 


N 


(i)/Z 

2 a 


N 




( 9 ) • 


Using (9) , (8) can be rewritten as 

E x (co) = (r s ( <»>) + j X s ( o>) ) H y ( g>) 


( 10 ) 


where R s is the surface resistance and X s is the surface reactance. 
Consider rewriting (10) as 


E x (g)) 


(r s (w) + jo>L s (G>))H y (6>) 


(ID 


with the resistance and inductance defined by 


R s (o» = 


6>/Lt 

nJ 2a 


( 12 ) 


N 


M«> 


2 a oi 



6 


To remove the frequency dependence of the surface resistance and 
inductance, these quantities are evaluated at a particular 
frequency and are subsequently treated as constants. Equation (11) 
then becomes 

E x (g>) = (r s + j °) L s ) H y ( o>) ( 13 ) • 


This is the required frequency domain constant surface impedance 
boundary condition. To incorporate this boundary condition into 
the FDTD algorithm, the time domain equivalent of (13) must be 
obtained. Performing an inverse Fourier transform operation on 
(13) results in 

E x (t) = R s H y (t) + L s A Hy (t) (14). 


This equation defines the time domain FDTD constant surface 
impedance boundary condition. 

To implement this constant surface impedance boundary 
condition, space and time are quantized by defining 

z - (k<Sz) - (k) (15) # 

t =* (n<St) =* (n) 


The Faraday -Maxwell law is then used to obtain the H y component in 
the free space cell next to the impedance boundary. Since the 
impedance boundary condition requires that the electric and 
magnetic fields are co— located in space and time, we assume that 
the magnetic field 1/2 cell in front of the impedance boundary and 
1/2 time step previous is an adequate approximation. The Faraday- 
Maxwell law yields 


-(/x 0 <Sx<Sz) 


d(nSt) 


HY n (k+l/2) 


= EX n (k + l) Sx - EX n (k) <Sx (16) 


Note the component EX n (k+l) of (16) is the electric field component 
at the impedance boundary. Quantizing space and time in (14) and 
using the result to eliminate EX n (k+l) in (16) gives 


-(m 0 Sz ) 


5 (n5t) 


HY n (k+l/2) 


R HY n (k + 1/2 ) + L 


5 (nfit) 


HY n (k + 1/ 2 ) 


- EX n (k) 
(17) • 


Notice that the HY n (k+1/2) term in (17) is time indexed at time 
step n. This term is approximated as 
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HY n (k + l/2) » - [HY n+V2 (k + l/2) + HY n 1/2 (k + 1/ 2) J (18). 

2 

Using (18), and approximating the time derivatives on the magnetic 
fields in (17) as finite differences gives 

-(ji 0 <S z + L s )( HY n+1/2 (k + 1/2 ) - HY n1/2 (k + l/2)) = — ( HY n+1/2 (k + 1/ 2 ) -HY n ~ 1/2 (k + 1/2 ) ) 

- 5tEX n (k) (19) . 

Solving for HY n+1/2 (k+1/2) in (19) yields 


HY n+1/2 (k + 1/2) 


/i n 5z +L, -R,5t/2] n 

= _2 ! ! HY n_1/2 (k + 1/2) 

/x 0 <5z +L S +R s 6t/2 


St 

H 0 Sz +L s + R s 5t/2 


EX n (k) 
( 20 ) . 


This equation implements the constant surface impedance boundary 
condition in the FDTD method. 

IV. FDTD Dispersive Surface Impedance Implementation 

To derive a similar relation to (20) valid over a wide 
frequency band, we begin with the same set of underlying 
assumptions as for the constant surface impedance. The primary 
exception is that the surface impedance will vary with frequency 
and will not be approximated by its value at a particular 
frequency. All frequency domain information is inversed Fourier 
transformed to equivalent time domain form. The SIBC is then 
implemented in the FDTD method with the required convolution using 
a recursive updating technique. 

The standard first order impedance boundary condition remains 
unchanged and is given by (8) . In a similar fashion as Tesche [9] , 
(8) is rewritten as 

z s (° ) )l 

E ( G>) = jo H ( o>) 

jw J y 


( 21 ) . 
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Defining 


/ 2 ( to) 

Z s (co) = — 

D&> 


( 22 ) 


and substituting (9) into (22) gives 


Z s (<*>) = 


(23) 


\] j ua 


Substituting (23) into (21), a modified surface impedance boundary 
condition is obtained as 


E ( (ji) = Z s (o>) [ jwH (o>) ] 


(24) . 


The time domain equivalent of (24) is obtained via an inverse 
Fourier transform operation as 


E x (t) 


Z„(t)» 




(25) 


where the asterisk denotes convolution, 

E x (t) = ^•• 1 [E x (o)) ] 

H y (t) = ^ 1 [ H y ( o>) ] (26) 

z'(t) = 


and the denotes the inverse Fourier transform operation. Note 

in (25) that as a -*■ OO , the boundary condition becomes E x (t)=0.0, 

which is required for a perfect conductor. To determine Z ? ' (t) , 
the Laplace transform variable s=jw is used in (23) to obtain 


Z s (s) = 


N y/s 


( 27 ) . 
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Using the Laplace transform pair [11] 


1 




(28) 


where the denotes the inverse Laplace transform operation; and 
Z s '(t) is then determined to be 


Zs(t) 


& A [ Zg ( o>) ] = \ not 


0, 


t >0 
t <0 


(29) . 


This is the required time domain surface impedance function. 
Substituting (29) into (25) and discretizing space and time 
according to (15) gives 


- HY n ( k + 1 / 2 ) 

3(n<St) 


EX n (k + l) 


\| no (n<St) 


Substituting (30) into (16) yields 


-M 0 5 z 


d(n£t) 


HY n (k+l/2) 


M 


7ra(n<St) 


3(n5t) 


■ HY n (k + 1/ 2 ) 


EX n (k) 

(31) 


The convolution in (31) is expressed as a summation to obtain 


-H Q Sz 


d (n<5t) 


HY n (k + 1/ 2 ) 



n-1 

= 

fJL <5t 


HQ m4) 


d ( (n-m) St) 


• HY n " fl1 (k + 1/ 2 ) 


Z 0 (m) - EX n (k) 
(32) 


where Z 0 (m) is the discrete impulse response. The discrete impulse 
response is obtained by assuming the fields are piecewise constant 
in time as 

m+1/2 

Z 0 (m) = f -Ida (33) 

m-i/2 ya 



10 


If m=0, the lower limit in (33) is 0. Approximating the time 
derivatives on the magnetic fields in (32) as finite differences 
results in 


HY n+1/2 (k + l/2 ) -HY n 1/2 (k + 1/2 ) j = -Z 1 ElHY n - fl ’ +1/2 (k + l/2) -HY nHn_1/2 (k + 1/2 ) 


Z 0 (m) 


m=0 


+ EX n (k) 




(34) 


where 


M 0 5z\| 


ust 

no 


(35) 


Equation (34) is suitable for computer implementation and includes 
the full convolution with all past field components. This full 
convolution would be impractical for large three dimensional 
problems; thus it is desirable to obtain a more efficient 
implementation. The development of a recursive implementation is 
the subject of the following section. 

V. Recursive Implementation 

Recently, Luebbers et. al. [3] extended the FDTD technique to 
dispersive dielectric materials using a time domain susceptibility 
function for polar dielectrics. In that paper, the time domain 
susceptibility function was a decaying exponential which permits 
the convolution summation to be recursively updated, thus avoiding 
the need for the complete time history of field components. Upon 
further examination of (34), it is clear that if Z 0 (m) can be 
approximated by a series of exponentials, then the SIBC can be 
efficiently evaluated using recursion. Figure 2 shows Z 0 (m) versus 
m, and it is clear that it can be approximated by a series of 
exponentials. Z 0 (m) is approximated as 

ZoW-Eay-" < 36) 

i =1 

where N is the number of terms in the approximation. One of the 
most accurate methods for obtaining an exponential approximation to 
an exact function or to a data set is Prony's method [10]. Figure 
2 also shows the Prony approximation to Z 0 (m) with N=10 and it is 
clear that N=10 provides an adequate approximation. Thus, using 
(36) with N=10 in (34) gives 


11 


10 n-1 r 

53 53 |HY n ' ro+1/2 (k + 1/2 ) - HY nHI,_1/2 (k + 1/ 2 ) 
EX n (k) 


HY n+V2 (k + 1/2) =HY n1/2 (k + 1/ 2 ) - 


a m 

a ; e 1 + 


1 + Zq ( 0 ) j = i m=1 

St 


where 


M 0 6z(1 + Z 1 Z 0 (O)) 
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Z 0 (O) = E a i 


(37) 


(38) 


i =1 


The convolution can now be recursively updated (see Appendix) to 
give 


10 


HY n+1/2 (k + 1/2) = HY n ' 1/2 (k + 1/2 ) 53tf(k+l/2) 

l+Z^Q (0) i= i 


(39) 


St 


ju 0 <Sz(l + Z 1 Z 0 (O) ) 


EX n (k) 


Note that only one past value of magnetic field is required to 
update the convolution summation. 

VI. One Dimensional Demonstration 

To demonstrate the constant and recursive FDTD SIBC, (20) and 
(39) were implemented in a one dimensional total field FDTD code 
for the geometry shown in Figure 1. The problem space size is 3 01 
cells, the impedance boundary is located at cell 3 00, and the 
electric field is sampled at cell 299. The maximum frequency of 
interest for each problem was 10 GHz. The incident electric field 
is a Gaussian pulse with maximum amplitude of 1000 V/m and has a 
total temporal width of 256 time steps. The frequency response of 
the incident pulse contains significant information to 12 GHz. Two 
computations were made with <7=2.0 S/m and a=20.0 S/m. The loss 
tangents at 10 GHz are 3.599, 35.99, respectively. The 

permittivity and permeability for the lossy dielectric were those 
of free space. The cell size and time step were 750 jum and 2.5 
psec , respectively. A tie point of 5.0 GHz was chosen for the FDTD 
constant SIBC. For each FDTD computation, a reflection coefficient 
versus frequency was obtained by first dividing the Fourier 
transform of the scattered field by the transform of the incident 
field at cell 299. The incident field was obtained by running the 
FDTD code with free space only and recording the electric field at 
cell 299. The scattered field is then obtained by subtracting the 
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time domain incident field from time domain total field. The 
results are compared with the analytic surface impedance reflection 
coefficient computed from 


z g (<»>) - n 0 

' z g (u>) + n 0 


(40) 


where Z s (<i)) is given by (9) and 'Ho is the free space wave 
impedance. The phase of the FDTD reflection coefficient was 
corrected to account for the round trip phase shift of one cell 
since the FDTD reflection coefficient is computed from electric 
fields recorded one cell in front of the impedance boundary. 


The high conductivity surface impedance of (9) is an 
approximation to the general surface impedance for lossy 
dielectrics given by 


Z s (o>) 


j<*>M 


\J a + jcoe 


(41) . 


The advantage of using (9) over (41) for the FDTD SIBC 
implementation is that the resulting time domain impulse response 
is independent of the conductivity. The exponential approximation 
needs to be performed only once and not each time the conductivity 
is changed. 

Figures 3-4 show the FDTD constant and recursive SIBC 
reflection coefficient magnitude and phase results versus the 
analytic SIBC results for a= 2.0 S/m. Notice the agreement between 
the curves is good, and the maximum error is about 0.02 at 10 GHz 
in Figure 3 . 

Figures 5-6 show the FDTD constant and recursive SIBC 
reflection coefficient magnitude and phase results versus the 
analytic SIBC results for a=20.0 S/m. Notice the agreement between 
the curves is excellent. 

Since the FDTD SIBC implementation is an approximation to an 
analytic SIBC, some amount of divergence between the SIBC curves 
and the analytic SIBC solution is to be expected with increasing 
frequency. As frequency increases, the effective number of cells 
per wavelength decreases and the FDTD SIBC becomes a rougher 
approximation to the analytic SIBC. To observe this error trend, 
the same one-dimensional test problems as above (using the 
dispersive SIBC only) were reevaluated with larger cell sizes equal 
to twice and four times the original cell size. This is equivalent 
to having 20 and 10 cells/X 0 in the free space region, 
respectively. Figures 7 and 8 show the FDTD dispersive SIBC 
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reflection coefficient magnitude and phase results versus the 
analytic SIBC for u= 2.0 S/m using the original cell size and the 
larger cell sizes. Notice that for each increase in cell size, the 
agreement between the SIBC curve and the exact solution is reduced 
by a factor of two. This indicates that the error in the SIBC 
implementation is 0(5z) over the range of cell sizes examined here. 
The constant SIBC exhibited similar agreement reductions at the 5 
GHz tie point for larger cell sizes. 

VIII. Two Dimensional Demonstration 

As a practical application of the FDTD dispersive SIBC, 
frequency domain scattering width was computed from an infinite 
square cylinder for two scattering angles, 0=0.0 and 0=30.0 degrees 
using a full two dimensional TM scattered field code. The cylinder 
was 0.99 cm square and had parameters e=e 0 , M=M 0 / and a=20.0 S/m. 
To illustrate the applicability of the SIBC, the cylinder was 
modeled in two ways. The first was a normal FDTD computation with 
a grid size of 10 cells/A (at 10 GHz) inside the conducting 
cylinder and the second was a SIBC computation with a grid size of 
10 cells/A. in free space (at 10 GHz) . Figure 7 shows the two 
dimensional field components and the cylinder dimensions (in cells) 
for the FDTD and SIBC computations. For the FDTD computation, the 
cylinder was modeled using 198 cells in the x and y directions, the 
cell size was 500 /im, and the time step was 1.18 ps. For the SIBC 
computation, the cylinder was modeled using 32 cells in the x and 
y directions, the cell size was 0.003 ra and the time step was 7.07 
ps. For both computations, a 100 cell border between the cylinder 
and the absorbing boundary was chosen, the total number of time 
steps was 1024, and an incident Gaussian pulse with total pulse 
width of 64 time steps was chosen. The near zone fields were 
transformed to far zone fields by a two-dimensional near zone to 
far zone transformation [12]. 

Figure 10 shows the scattering width magnitude versus 
frequency for a scattering angle of 0=0.0 degrees using the FDTD 
computation and the SIBC computation. Notice the good agreement 
over the entire frequency bandwidth for the dispersive SIBC. 

Figure 11 shows the scattering width magnitude versus 
frequency for a scattering angle of 0=30.0 degrees using the FDTD 
computation and the SIBC computation. Notice again the good 
agreement over the entire frequency bandwidth for the dispersive 
SIBC. 


IX. Summary 

One dimensional FDTD implementations of constant and 
dispersive surface impedance boundary conditions have been 
presented. The corresponding time domain impedance boundary 
conditions have been derived and their validity demonstrated by 
one-dimensional computation of the reflection coefficient at an 
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air-lossy dielectric interface at a single frequency and over a 
wide frequency bandwidth. The applicability of the SIBC to two- 
dimensional scattering problems was demonstrated by scattering 
width computation from an infinite square cylinder. For both the 
one and two dimensional cases, the dispersive FDTD results were 
shown to be in good agreement with exact results over the entire 
bandwidth. Considerable computational savings were illustrated and 
a recursive updating scheme was implemented which permits efficient 
application of a dispersive surface impedance boundary condition to 
practical scattering problems. 

Future extensions of this surface impedance concept currently 
under investigation are implementation in three dimensions, 
inclusion of surface curvature, dispersive dielectric and magnetic 
materials and thin material layers. 
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Appendix 

The purpose of this Appendix is to show how the discrete 
convolution of the SIBC in (37) can be done recursively. The 
convolution in (37) is 

1° n-1 r 1 

53£[HY nH " +1/2 (k + l/2) -HY n " n ' 1/2 (k + l/2)j a 1 e a,m ( 42 ) 

i=1 m=1 


Consider n=2 , and (42) becomes 


10 / \ 

£\HY 3/2 (k + l/2) - HY 1/2 (k + l/2))a j e a ' 


i =1 


(43) 


Now define 


(k + 1/2) 


( HY 3/2 (k + 1/ 2 ) - HY 1/2 (k + l/2))a,.e a ' (44) 


Next for n=3, (42) becomes 


10 2 


52 53 (HY 7/2 (k + l/2 ) - HY 5/2_m (k + l/2) ) a^ ' 


i =1 m=1 


(45) 


Expanding (45) gives 


10 / 

E(> 


HY 5/2 (k + 1/2 ) - HY J, ‘(k 


3/2, 


: + l/2))a i 


a 

e 1 + 


i=i 

(hy 3/2 (k + 1/2 ) - HY 1/2 (k + 1/2 )) aje 24 ' 


(46) 


Substituting (44) into (46) we obtain 
10 


52 ( HY S/2 (k+1/2 ) - HY 3/2 (k+1/2 jlaje 0 ’ + e"’^ (k + 1/2 ) 


i =1 


Equation (47) can be generalized for any time step n as 

10 17 \ n-1 

52 (hy 0 * 172 (k + 1/2 ) - HY n_3/2 (k + 1/2)) a-e*' + e a, i|;" (k + 1/2) 


i =1 


(47) 


(48) 


with 
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Vi* (k+1/2) 


(hy"~ 1/2 (k + 1/2) -HY n ' 3/2 (k + 1/2)) 


a,e 


e a 'i|i n ‘ 1 (k + 1/2) 


(49) 


and 


(k + 1/2) = i|»° (k+1/2) = 0.0 


(50) 
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FIGURE TITLES 

1. Problem geometry showing one— dimensional FDTD grid and planar 
free space-conductor interface. 

2. FDTD dispersive SIBC discrete impulse response Z Q (m) versus m 
and Prony approximation using 10 terms. 

3. Reflection coefficient magnitude versus frequency for normal 

incidence plane wave calculated for a=2 . 0 S/m using FDTD 

constant and dispersive SIBC and analytic solution. 

4. Reflection coefficient phase versus frequency for normal 

incidence plane wave calculated for a=2 . 0 S/m using FDTD 

constant and dispersive SIBC and analytic solution. 

5. Reflection coefficient magnitude versus frequency for normal 
incidence plane wave calculated for a=20.0 S/m using FDTD 
constant and dispersive SIBC and analytic solution. 

6. Reflection coefficient phase versus frequency for normal 

incidence plane wave calculated for ct= 20.0 S/m using FDTD 
constant and dispersive SIBC and analytic solution. 

7. Reflection coefficient magnitude versus frequency for normal 

incidence plane wave calculated for a=2 . 0 S/m using FDTD 

dispersive SIBC with original and larger cell size and 

analytic solution. 

8. Reflection coefficient phase versus frequency for normal 

incidence plane wave calculated for a=2.0 S/m using FDTD 

dispersive SIBC with original and larger cell size and 

analytic solution. 

9. Two dimensional geometry for scattering width computations 
from an infinite square cylinder with a=20.0 S/m using normal 
FDTD and FDTD dispersive SIBC. 

10. Scattering width magnitude versus frequency at scattering 
angle <£=0.0 degrees from an infinite square cylinder with 
a=20.0 S/m using normal FDTD and FDTD dispersive SIBC. 

11. Scattering width magnitude versus frequency at scattering 
angle 0=30.0 degrees from an infinite square cylinder with 
a=20.0 S/m using normal FDTD and FDTD dispersive SIBC. 



FDTD SIBC Discrete Impulse Response Zo( 
Exact and Prony approximation 



{(m)°z } 3oi 


100 150 200 250 300 350 400 450 500 




ient magnitude versus frequency 



a jo sprniuSew 


Frequency (GHz) 





Reflection coefficient phase versus frequency 
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Reflection coefficient phase versus frequency 
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Scattering width magnitude versus freq 
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Scattering width magnitude versus frequency 
TM cylinder, 0.99 cm square, a=20.0 S/m, < 2 ? = 30.0 degrees 



(ra/gp) apnq.mSiera qypui Suijayyeog 



Frequency (GHz) 



